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■ Abstract: We review some recent results concerning integrable quantum field theories in 

[ — ■ 1 + 1 space-time dimensions which contain unstable particles in their spectrum. Recalling 

first the main features of analytic scattering theories associated to integrable models, 
we subsequently propose a new bootstrap principle which allows for the construction of 
particle spectra involving unstable as well as stable particles. We describe the general 
' Lie algebraic structure which underlies theories with unstable particles and formulate 

a decoupling rule, which predicts the renormalization group flow in dependence of the 
relative ordering of the resonance parameters. We extend these ideas to theories with 
an infinite spectrum of unstable particles. We provide new expressions for the scattering 
amplitudes in the soliton-antisoliton sector of the elliptic sine-Gordon model in terms of 
infinite products of q-deformed gamma functions. When relaxing the usual restriction 



Oh, on the coupling constants, the model contains additional bound states which admit an 

interpretation as breathers. For that situation we compute the complete S-matrix of all 
sectors. We carry out various reductions of the model, one of them leading to a new type 



of theory, namely an elliptic version of the minimal SO(n)-affine Toda field theory. 



1. Introduction 

The structure of integrable quantum field theories (IQFT) in 1+1 space-time dimensions 
has been unravelled to a very large extend. Many theories can be solved even exactly, 
that is to all orders in perturbation theory, in this context. However, the large majority of 
investigations concentrates on theories which involve exclusively stable particles, despite 
the fact that in nature most particles are unstable. Since of course one of the motivations 
to study IQFT is to reproduce realistic features, there is an apparent need to investigate 
also theories which have unstable particles in their spectrum. The aim of this talk is to 
review some recent results which deal with such theories. 



°Proceeding of the workshop on "Infinite dimensional algebras and quantum integrable systems" (Faro, 
Portugal, July, 2003). We thank the organizers, especially Nenad Manojlovic, for their kind hospitality and 
untiring engagement to make things work. 
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2. Analytic scattering theory of factorizable integrable models 

Since not all participants of this conference work directly on integrable quantum field 
theories, we briefly recall some well known facts on analytic scattering theories in 1+1 
space-time dimensions. Having in mind to emphasize features related to unstable particles 
this will also be useful to the experts. As a starting point in every scattering theory one 
requires a complete set of asymptotic in and out states (t — ► ±00). These states consist of 
operators Z^{p) acting on the vacuum |0) and creating in this way a stable particle of the 
type \i with momentum p. Already at this point enters the fundamental difference between 
stable and unstable particles. Even though experimentally unstable particles with a very 
long lifetime can very often not be distinguished from stable ones, mathematically they are 
very distinct. They can never be associated to an asymptotic state, even when they have 
an extremely long lifetime, as by their very nature they will have decayed in the infinite 
future or were never produced in the infinite past. Then the scattering matrix is defined 
to be the operator which relates a stable n-particle in state to a stable m-particle out state 

z,je' m ) . . . z^Oi) |o) out = s^:;^, . . . e n )z^e x ) . . . z^n) |o) in . (2.1) 

Conveniently one parameterizes the two-momentum by the rapidity 8asp = m(cosh 9, sinh 9) 
Now there are some very special features happening in integrable (that means here there 
exists at least one non-trivial conserved charge) quantum field theories in 1+1 dimensions 
HI |2|, [3|, [|, H . There is no particle production and furthermore the incoming and outgoing 
momenta coincide 

{e' 1 ,e' 2 ,...e' m } = {0 1 ,e 2 ,...e n } with n = m . (2.2) 

In addition, the n-particle S-matrix factorizes into a set of two-particle S-matrices 

^&:&(C--*») = n W'*^) • ( 2 - 3 ) 

l<i<j<n 

Obviously, this is a considerable simplification in comparison with the general situation 
( |2.1| ), as it means that once we know the two-particle 5-matrix, we control the entire 
scattering matrix. Because of this fact, we refer from now on to the two-particle scattering 
matrix as the .S-matrix. 

How does one construct this S-matrix? In general one is limited to the use of perturba- 
tion theory in the coupling constant. In particular in higher dimensions that is essentially 
the only method available. In contrast, two dimensions are very special as they miracu- 
lously allow to determine S exactly to all orders in perturbation theory. This is one of the 
major successes of this area of research and one of the reasons for the continued interest in 
such theories. The original ideas which lead to explicit expressions for S go back to what 
is called the bootstrap approach ||, ||, §]. It consists of using various properties for the 
scattering matrix, which one motivates by some physical principles in order to set up an 
axiomatic system for S in the hope that it might be so constraining that it determines S 
completely. Indeed, these hopes are not in vain. 
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We recall the S'-matrix properties: 

i) Lorentz invariance 

Dealing with relativistic scattering theories, we expect the scattering matrix to be 
Lorentz invariant, i.e. it should depend only on covariant combinations of the momenta. 
The Mandelstam variables are precisely such quantities, see e.g. 0. In 1+1 dimensions 
only one of them is independent, usually taken to be s a b = (p a + Pb) 2 = m„ + m\ + 
2m a m b cosh(0 ' a — 9 b ). Hence, Lorentz invariance is simply guaranteed when S depends 
either only on s ab or the rapidity difference 9 ab := 9 a — 6 b 

Sab(p a ,Pb) = S ab (9 a ,e b ) = S ab (s ab ) = S ab (9 ab ). (2.4) 

Since s ab admits the interpretation as the total energy in the centre of mass system, 9 ab 
has to be real for a physical process, such that s ab > (m a + m b ) 2 . 

ii) Hermitian analyticity 

As a central assumption of analytic S'-matrix theory Q one assumes that the S-matrix 
can be continued to the complex plane and depends on s ab ,9 ab £ C. Physical scattering 
amplitudes are then assumed to be real boundary values of analytic functions, which can 
be obtained from a generalization of Feynman's ie prescription of perturbation theory 

physical = ^ ^ + i£) = Sab{9) seR ,e,9e R+. (2.5) 

The choice of the signs is important and relates to causality. Since a two-particle wave- 
function, having here plane waves in mind modulated by some enveloping function, will 
depend on the sum of the momenta, i.e. on y/s ab , one has lost the single valuedness of the 
scattering matrix by an analytic continuation. This is remedied by branch cuts along the 
real axis at s ab > {m a + m b ) 2 and s ab < (m a — m b ) 2 , the latter being motivated by crossing 
see iv). Hermitian analyticity is now a postulate which states how to continue over these 
cuts g I 

lim S ab (s + ie) = lim S ab (s - ie) & S ab (9) = [S ba (-9*)}* . (2.6) 

e^O e— >0 

once more for s£l,£,S£ R + . The equivalence is due to the fact that the analytic con- 
tinuation s + ie «-> s — ie corresponds to 9 <-> —9. Often one merely uses real analyticity 
S a b(9) = [S a b(— #*)]* instead of (|2.6| ), which only coincides when S ab = S ba , that is in 
parity invariant theories. This difference is very important with regard to the theories con- 
sider below, which involve unstable particles as they unavoidably break parity invariance. 



Further support for (2J3) comes from perturbation theory ||], general considerations in 
analytic S-matrix theory |l~0| , ^] and explicitly constructed examples, 
iii) Unitarity 



Assuming that the states in ( |2.1| ) are complete and orthogonal, the operator which 
maps them to each other has to be unitary 

SS^ = S^S = 1 . (2.7) 



The combination of (2.6) and (2/7) leads to the simpler relation S ab (9)S ba (—9) = 1, 
which may also be derived from applying twice the Zamolodchikov algebra Z a {9\)Z b {92) = 
S ab{8 12) Z b (9 2) Z a {9\). 
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iv) Crossing symmetry 

Crossing symmetry can be motivated by the Lehmann-Symanzik-Zimmermann (LSZ) 



formalism [11] and consists of the replacement of an incoming particle a by its anti-particle 



a with reversed momentum. A discussion of the anti-particle theorem can be found in 



|1C[. The prescription amounts to the continuation of the Mandelstam variable s ab to the 
variable t ab = (p a - p b ) 2 

lim S ab {s + ie) = lim S b - a {t - ie) ^ S b - a {9) = S ab (iir - 9) . (2.8) 

It is easy to check that the analytic continuation s + ie «-> t — ie corresponds to 9 «-> in — 9. 
v) Yang-Baxter equation 



In (|2.3| ) we already indicated that the conserved charge(s) of an integrable theory can 
be used to disentangle an n-particle scattering process into a consecutive scattering of two 
particles only. An additional consequence of this argument is that the order in which this 
takes place does not matter, such that two different orderings are taken to be equivalent. 
As in general the S-matrices do not commute, this leads to a new constraint. In other 



words this amounts to say that the operators Z in (2.1) obey an associative algebra. As a 



result of this one obtains the Yang-Baxter equation [12, [I 



S(e 12 ) ® 5(0i3) ® S(6 23 ) = S(9 23 ) ® 5(0i3) ® 5(012) • (2.9) 

For diagonal theories, i.e. when backscattering is absent, we simply have S^f(9) — ► S ab (9) 
such that fl2.9|) is trivially satisfied, 
vi) Fusing bootstrap equation 

By the same reasoning as in v) integrability, i.e. factorizability, yields a further con- 
straining equation, when two particles are allowed (what that means is discussed in vii)) to 
form bound state. For instance, the particles a, b fuse to a third particle c, i.e. a + b — > c. 
One makes now a further assumption, sometimes referred to as nuclear democracy, namely 
that also the particle of type c exists asymptotically. Then, by integrability, it is equiva- 
lence if a third particle, say I, scatters with the bound state c or consecutively with the 
two particles a, b. For 5 this reads 

Sic(0) = S la (9 + iv b ac )S lb (9 - if)l) . (2.10) 

The fj b ac S R + are the fusing angles specific to the individual theory considered. It is clear 
that the assumption of nuclear democracy does not hold if c is an unstable particle, such 
that (ElcD can not be valid in the form stated for that case. We will now indicate the 



origin for the possibility to form bound states, which is the 
vii) Pole structure 

In general, the 5-matrix can have a quite intricate singularity structure consisting of 
poles of finite order distributed all over the complex s, 0-plane. A strong further constraint 
is to assume that all singularities which emerge in 5 admit a consistent explanation. As 
a slightly weaker assumption one could suppose that all explainable poles form a coherent 



system, in the sense that the bootstrap (2.10) closes etc., and allow some redundant poles. 
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Single order poles are most important as they determine the particle spectrum of the 
theory. In the s-plane they might be on the real axis between the two branch cuts at 
s = m§, interpreted as an on-shell bound state particle, or in the second Riemann sheet 
at s = (nic — iT s /2) 2 corresponding to an unstable particle with finite lifetime r = l/r g . 
The discussion is more conveniently carried out in the 0-plane, since S(6) is a meromorphic 
function unlike S(s). Near the singularity S has to be of the form 

S ab (9) ~ lRkb . (2.11) 
(9 - trfa + a c ab ) 

Depending on the location and signs of the residues we have the following interpretations 

s-channel bound state: R c ab G R + ,rj ] b G R + ,a c ab = 
t-channel bound state: R c ab G R~,rj a b G R + , a c ab = 
unstable particle: R c ab G R, ff ab G R~,cr c ab G M~ 

The relation between the poles in the s and 8 planes are the Breit-Wigner (BW) equations 



Ml 



r 2 

m 2 j- = m 2 + mi + 2m a m\ ) cosh a c b cos rf ab (2-12) 

c 4 

m B T s = 2m a mb sinh a c ab sin ?7^ b , (2-13) 

which allow to express the mass m_ and decay width V- of the unstable particle as functions 
of m a , ma, 77^, For the stable particle formation we have the following relation between 
the fusing angles in ( p. 10 ) and the poles in ( 2.11| ): fj b ac = ir — rf ab . Note further that for the 



unstable particle formation in ( [2.11 ) we made the definite choice that the unstable particle 
c is formed in the process 

a + b^c (2.14) 

rather than in the not equivalent one b + a. It is clear that parity has to be broken, as with 
the choice rf ab ,o- c ab G R~ the amplitude Sj,a(G) will have a pole at iir — irj c ab + a c ab , leaving 



in ( 2.13 ) the choice that either m s < or T E < 0, which is of course both non-physical. 

Below we will be particularly interested in the situation for large resonance parameters 
a c ab , when the mass of the unstable particles can be approximated as 

mc ~ ^/m a mbe~ CTab//2 . (2-15) 

In terms of perturbation theory in a coupling constant /3, the pole ( |2.11| ) would be of 
second order, i.e. R^p 2 ), corresponding to a tree diagram. Similarly, higher order poles 
admit interpretations in form of more complicated singular Feynman diagrams. In some 
simple theories, such as for example sine-Gordon, the highest order of the poles is two. 
In that context it was suggested [15] that such type of poles are of order ft 4 and may be 
viewed as box diagrams. For quite some time higher order poles were ignored and also 
here we will not enter into a deeper discussion of them, which can be found for instance in 



1 16, 17]. From what is said it is clear that such poles will not alter the particle spectrum. 



Nonetheless, one should be able to draw the relevant Feynman diagrams, which means one 
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needs certain three-point couplings to be non-vanishing. Consequently this is a constraint 
on the existence of certain three point couplings R c ab . 

Remarkably, the constraints i)-vii) allow to determine the S -matrix exactly, that is to 
all orders in perturbation theory. However, one should say that the solution constructed 
this way is not unique, as there exists always the possibility to multiply with so-called 
CDD-factors [18]. To fix them requires additional arguments beyond the scheme outlined 
above, such as ultraviolet limits, certain inputs from Lagrangians, etc. 

2.1 A proposal for a construction principle of unstable particle spectra 

We have seen in the previous section, that there exists a powerful construction principle for 
the spectrum of stable particles, consisting of solving the equations (axioms) i)-vii). For 
unstable particles we do not have yet such a construction tool, as by now they emerge rather 
passively as poles in the unphysical sheet as by-products in the scattering process of two 
stable particles. Furthermore, a description of the scattering process of an unstable particle 
with another stable or unstable particle is entirely missing in this context. Obviously, 
scattering processes involving unstable particles do occur in nature, such that the quest for 
a proper prescription is of physical relevance. In addition, one aims of course always at a 
description which has predictive power. 

From what has been said, it is clear that such a description can not be a scatter- 
ing theory in the usual sense, since for that one requires the particles involved to exist 
asymptotically. Any unstable particle will vanish in this limit rendering such formulation 
meaningless at first sight. Nonetheless, some particles have extremely long lifetimes, and 
seem to exist quasi infinitely long from an experimentalists point of view. It appears there- 
fore natural to seek a principle closely related to the conventional bootstrap for stable 
particles. Inspired by this we proposed the following construction principle: 

Let us assume that in the time interval < t < t s we can formally associate to the 
unstable particle some creation operator z\(Q), with ]im t -^oa z\{0) = 1 if t 5 < oo. It is 
clear that these operators do not exist asymptotically, but for the stated time interval they 
can mimic an asymptotic state. Let us now further suppose that these operators satisfy a 
Zamolodchikov algebra 

z a (e 1 )z b (e 2 ) = s ab (e 12 )z b (6 2 )z a (9 1 ) (2.16) 
z a (Q x )z h (e 2 ) = s a b{0i2)z b (e 2 )z a {e{) (2.17) 

which can be used to generate an S-matrix type of amplitude S ab , describing the scattering 
of one unstable particle with a stable one ( [2.16| ) or the scattering of two unstable particles 



( 2.17 ). We may proceed as before and ask which type of properties might be satisfied for 
S. 

i) Lorentz invariance 



As already indicated in ( [2.16 ), ( |2.17 ) it is natural to expect Lorentz invariance also for 
this amplitude such that S depends only the rapidity difference 9 ab 

S a b{Pa,Pb) = S ab (9 a ,e b ) = S ab (s ab ) = S ab (6 ab ). (2.18) 
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ii,iii) Hermitian analyticity, unitarity 

We will not make any assumption on hermitian analyticity here and in fact we do not 
expect unitarity to hold, since the states formed with the Z are not complete. However, 
applying ( 2.16| ) or ( 2,17| ) twice yields 



Sab(0)§ ba (-0) = 1, 



(2.19) 



which also holds for S, derivable from combining (^1|) and (^^) in that case. In fact, also 
for the construction of S it is really only the corresponding equation to ( |2.19| ) which is 
employed, rather than individually ( |2.6| ) and ( |2.7| ). 
iv) Crossing symmetry 

The validity of crossing can also be argued as before, but now we have to continue as 

s b - a (-e) = s ab (i7T + e) , (2.20) 



lim S a b{s 



is) 



lim Sba(t + ie) 



which in the #-plane amounts to the same equation as the one for S. 
v) Yang-Baxter equation 

Supposing the algebra related to ( 2.16| ), ( 2.17 ) is associative we have by the same 
reasoning as for stable particles the Yang-Baxter equation 



S(0 2 3) ® 5(013) ® S(0i 2 ) 



(2.21) 



S(9u) ®S(6 13 )®S(6 23 ) 

vi) Fusing bootstrap equation 

We commence with the fusing of two stable particles to create an unstable particle as 
in the process (EQ3) . To this process we can associate bootstrap equations almost in the 
usual way. We scatter for this with an additional stable or unstable particle, say of type I, 
and obtain the S bootstrap equations 



s la (e-f ca )s lb (e + -fi) =s rc (B), 



(2.22) 



where 7 = ±in— 7, 7 = irj — a and also 7 — > —7 is not a symmetry. The angles should 
be measured anti-clockwise, which explains the signs. We also note that we do not assume 
parity invariance, such that in general 7^ a 7^ 7^. With the help of ( p. 19 ), (2.20) one derives 
the bootstrap equations for the opposite parity and the ones for the crossed processes 
a + c — ► b and b + c — > a and from ( [2.22 ) 

Sci(d) = Sai(0 + tL) Su(e - to) » (2-23) 

Su(e) = sue - to) sue ±™- t a - t c ), (2.24) 
Sim = sue + ta)~sue ±™ + 7 L + to) ■ (2.25) 



From the crossing relation for the "scattering matrix" and (2.24) or ( p. 25 ) one obtains 
some relations between the various fusing angles 



tb + ta + to 



±iir . 



(2.26) 



At first sight this looks very much like the usual bootstrap prescription, but there are 
some differences. As is clear from the scattering process of two stable particles producing 
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an unstable one, the angle ^ c ah is not purely complex any longer as it is for the situation 
when exclusively stable particles scatter. As a consequence, this property then extends to 
the other angles 7^ and 7?, in (2.22), which also possess some non-vanishing real parts. 



Note that ( |2.26| ) implies that the real parts of the three angels involved add up to zero. 
At this point we do not have an entirely compelling reason for demanding that, but this 
formulation will turn out to work well. 

Of course the above equations are only a proposal, which needs to be put on more solid 
ground. Nonetheless, at this point our proposal gains support from self-consistency and 
its predictive power, which may be double checked: a) The bootstrap closes consistently 
for many non-trivial examples, which we calculated. As for stable particles this is never 
guaranteed and by no means self-evident, b) The bootstrap yields the amount of unstable 
particles together with their mass. This prediction can be used to explain a mass degeneracy 
of some unstable particles which can not be seen in a thermodynamic Bethe ansatz (TBA) 
analysis for the concrete example of the homogeneous sine-Gordon (HSG) models, see 
below, c) The bootstrap is in agreement with a general Lie algebraic decoupling rule, which 
we also present below, describing the behaviour when certain resonance parameters tend 
to infinity, d) The bootstrap yields the three-point couplings of all possible interactions, 
that is, involving stable as well as unstable particles. 

2.2 An example: The g^-HSG model 



The gfc-homogeneous sine-Gordon models (HSG) (2^, 21], with g being a simple Lie algebra 
of rank t and level k, will be our standard example in what follows. In fact they have been 
the first models with a well defined Lagrangian containing unstable particles which have 



been the subject of a systematic analysis |22|, |23|, |2j, |25|, |26|, 27, 28, 19, |29f| . They can be 
viewed as perturbed conformal field theories (CFTs) 1 

W Gfc -HSG = T-(-G k /U(iy-CFT ~ A / d2x( t>{x,t) . (2.27) 
The underlying ultraviolet CFT is a Wess-Zumino-Novikov- Witten-Gfc/C/(l)^-coset theory 



1 31, [3^]. The correspondig Virasoro central charge c is computed with standard argu- 
ments of |3^] and the perturbing operator eft is identified with a primary field of conformal 
dimensions A, A. One finds 

n kh-h y A A ^ 

c = £- — and A = A = - — . 2.28 

k + h y k + h v K ' 

Here (/i v ) h is the (dual) Coxeter number of g. For simplicity we will drop in the following 
the explicit mentioning of the subalgebra U(lY which were indicated in ( |2.27| ). 

The scattering matrix for g^-HSG-models with g simply laced algebras was constructed 



in [22]. For k = 2 it can be brought into the simple form 

Syfraij) = ^^(a,^)^ , 1 < i,j < I (2.29) 



1 For the particular case of the SU(3)2-HSG model it was shown ^5(J that it can be described alternatively 
as a perturbation of a tensor product of two minimal CFTs. 
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where I denotes the incidence matrix of g and e{x) is the step-function, i.e. e(x) = 1 for 
x > 0, e{x) = — 1 for x < 0. It is convenient to use the abbreviation 

(a,x) :=tanh(fl + f7-nrx/4)/2. (2.30) 

Let us now consider the concrete case SU (3)2- We can start with the known part 
of the scattering matrix ( 2.29| ) for the stable particles, and leave the remaining entries 



which involve unstable particles unknown. From this we construct consistent solutions 
to the bootstrap equations ( ^22|) , (|2~2l ) and d2~25l) . We can fix the imaginary parts of 



the fusing angles by the requirement that for vanishing resonance parameters we want to 
reproduce the masses predicted by the Breit-Wigner formula. Choosing the masses of the 
stable particles to be mi = 771,2 = m the one for the unstable results to tt7( 12 ) = \f2m. 
This argument does not constrain the real parts of the fusing angles, such that they are 
not completely fixed and still contain a certain ambiguity. The different choices of these 
parameters give rise to slightly different theories. First we consider the case CT21 > 0. 



Ct\ CK 2 

For this choice of the resonance parameter, we then find the following bootstrap equations 

%2)(0) = S n (e + (1 - u)a 12 + in/±)S l2 {6 - ua l2 - *7r/4) (2.31) 
from which we construct 

/ _i -(ai2,2) -((1- j/)*7 la ,3)\ 

S S U(3)(0,a 12 )= ((7 2 i,2) -1 -(^21,1) . (2.32) 
\-((v - l)a 12 ,l) -(i/<7i3,3) -1 J 

Here we label the rows and columns in the order {1,2, (12)}. According to the principles 



outlined above, the S-matrix ( |2.32| ) allows for the processes 

l + 2->(12), 2 + (12)->1, (12) + l-»2. (2.33) 
The related fusing angles are read off from ( |2.32| ) as 

(12) 2 ^* 7r 1 3i7r 

7l2 = -y + 7(12)1 = 4 - + l 1 - ^) fJ 12, 72(12) = 4- + V °U- V 2 - 34 ) 



and are interrelated through equation ( 2.26| ), which still holds even though the 7's have 



non- vanishing real parts. We can employ these fusing angles and compute the masses and 
decay widths by means of the Breit-Wigner formulae ( |2.12 ) and ( |2.13| ). Taking again for 



simplicity m\ = m 2 = m and in addition v = 1/2, we obtain for the first process in ( p.33| ) 

777(i2) = V^ttt- cosh o"2i /2 and T(i 2 ) = 2\/2tt7 sinh o"2i/2 . (2.35) 

Employing now also in the process 2 + (12) — * 1 the Breit-Wigner formula, we reproduce 
in the limit CT12 — * the values mi = m and Ti = 0. Likewise, in the last process in ( 2.33j ) 
we obtain m2 = m and T2 = 0. 
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The asymptotic limit t — > oo becomes meaningful when we operate on an energy scale 
at which the unstable particle has not even been created yet, i.e. rV 12 ) -^00 = cr 21 — ► 00. 
In that case the theory decouples into two SU(2)2-models, i.e. free Fermions, with Sn = 
5 2 2 = — 1- This is a simple version of the decoupling rule ( |3.3| ). 

Next we consider a different theory with a 12 > 0. 



a\ a 2 

Taking also in this case for simplicity v = 1/2, we find the following bootstrap satisfied 

Sj(i2)(0) = S l2 {6 - ai 2 /2 + m/A)S lx {0 + a 12 /2 - ztt/4), (2.36) 
which yields the S-matrix 

/ -1 (ai 2 ,2) -(a 12 /2,l)\ 

S S U(3)(0,O- 21 )= -((721,2) -1 -((721/2,3) . (2.37) 

\-((7 2 l/2, 3) -((71 2 /2,1) -1 J 

The S-matrix (|2.37j ) allows for the processes 

2 + 1 -> (12), 1 + (12) -> 2, (12) + 2-»l, (2.38) 



instead of (2.33). Now the fusing angles are read off as 



(12) ™ . 2 3i7F 0*12 1 3i7T CT12 ,„ „ m 

72i = - y + 7i(i2) = — 4 2", 72(12) = — 4^ 2" ^ ' 



and also satisfy ( 2.26| ). The masses and decay width are obtained again from ( [2,12 ) and 



( 2.13 ) with (7i2 — > o~2i- As a whole, we can think of this theory simply as being obtained 
from the Z2-Dynkin diagram automorphism which exchanges the roles of the particles 1 and 
2. However, since parity invariance is now broken this is not a symmetry any more and the 
two theories are different. In the asymptotic limit a~\ 2 — > 00, we obtain once again a simple 
version of the decoupling rule ( |3.3| ) and the theory decouples into two SU(2) 2 -models. 

The next example, SU(4)2-HSG, is more intriguing as it leads to the prediction a 
new unstable particle. Proceeding in the way as before we construct the corresponding 
amplitudes S, for details sec |l9| . We found there the processes 



1 + 2^ (12), (12) + 1- 2, 2+(12)-> 1, 

3 + 2^ (23), (23) +3 ^ 2, 2 + (23) -> 3, K ' 1 

which simply correspond to two copies of SU(3)2-HSG. It is interesting to note that the 
amplitudes S^ 2 )3 and £(23)1 contain poles at 

(123) _ 0-21 ~ 2(723 _ 3i7T (123) _ ^23 ~ gggl _ 3ivr 

'(12)3 ~ 2 4 7 23 ) 1 ~ 2 4 ' 



which yield the possible processes 



(12) + 3 -» (123), (123) + (12) 3, 3 + (123) -> (12), 
(23) + 1 -> (123), (123) + (23) -* 1, 1 + (123) -> (23). 



(2.42) 
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An interesting prediction results from the consideration of the first two processes in ( 2.42| ) 



Making in the first process the particle (12) and in the second the particle (23) stable, by 
02 — ► o\ and by 02 — ► 03, respectively, both predict the mass of the particle (123) as 

"1(123) ~ me kl3l/2 • (2.43) 

This value is precisely the one we expect from the approximation in the Breit-Wigner 
formula ( 2.1 5| ) . Note that in one case we obtain CT13 and in the other 031 as a resonance 



parameter. The difference results from the fact that according to the processes ( [2.42 ), the 



particle (123) is either formed as (1 + 2) + 3 or 3 + (2 + 1). Thus the different parity shows 
up in this process, but this has no effect on the values for the mass. 

In |l9| we presented more examples and remarkably we found consistency in each 
case. We take the closure of the bootstrap equations as a non-trivial confirmation for our 
proposal. 

3. Lie algebraic structure for theories with unstable particles 

There exist some concrete Lagrangian formulations for integrable theories with unstable 
particles, such as the aforementioned HSG-models ( 2.27] ). Inspired by the structure of 



these models, we present here a slighly more general Lie algebraic picture. We keep the 
discussion here abstract and supply below concrete examples. For our formulation we need 
an arbitrary simply laced Lie algebra g (possibly with a subalgebra h) with rank I together 
with its associated Dynkin diagram (see for instance [|33|]). To each node we attach a simply 
laced Lie algebra g« with rank li and to each link between the nodes i and j a resonance 
parameter aij = &i — o~j, as depicted in the following g/ h-coset Dynkin diagram 

&ij ®jk / 0~[ m 0~ mn 

• < * > * * / < * > •••• 



gl gi gj gfc Si / Si Sm gn 

Besides the usual rules for Dynkin diagrams, we adopt here the convention that we add an 
arrow to the link, which manifests the parity breaking and allows to identify the signs of 
the resonance parameters. An arrow pointing from the node i to j simply indicates that 
Uij > 0. Since we are dealing exclusively with simply laced Lie algebras, this should not 
lead to confusion. To each simple root of the algebras gj, we associate now a stable particle 
and to each positive non-simple root of g an unstable particle, such that 

# of stable particles = / # of unstable particles = . (3.1) 

i=i 

From the discussion above, we expect that the ex's will be associated to unstable particles, 
but we note that the 

£(£- 1) 

# of resonance parameters = (3-2) 

only agrees with the amount of unstable particles for h = £ + 1, e.g. for g = SU(£ + 1). 
Since the resonance parameters govern the mass of the unstable particles, this discrepancy 
is interpreted as an unavoidable mass degeneracy. 
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Concrete examples for this formulations are the g^-homogeneous sine-Gordon models 
|2C, 21], for which one can choose g to be simply laced and gi = • • • = g| = SU(k). This 
is generalized [34] when taking instead g to be non-simply laced and gj = SU(2k/af), 
with a.i being the simple roots of g. The choice gi = ... = g^ = g with g being any 
arbitrary simply laced Lie algebra gives the g|g-theories [35]. An example for a theory 



associated to a coset is the roaming sinh-Gordon model |36|], which can be thought of as 
g/ h = linifc^oo SU(k + l)/SU(k) with gi = . . . = g| = SU(2). It is clear that the 
examples presented here do not exhaust yet all possible combinations and the structure 
mentioned above allows for more combinations of algebras, which are not yet explored. 
One is also not limited to Dynkin diagrams and may consider more general graphs which 
have multiple links, i.e. resonance parameters, between various nodes. Examples for such 



theories were proposed and studied in [37|. 
3.1 Decoupling Rule 

Of special interest is to investigate the behaviour of previously defined systems when certain 
resonance parameters a become very large or tend to infinity. The physical motivation for 
that is to describe a renormalization group (RG) flow, which we shall discuss in more detail 
below. Here we present first the mathematical set up. 

Decoupling rule: Call the overall Dynkin diagram C and denote the associated Lie 
group and Lie algebra by Gq o- n d gc, respectively. Let Oij be some resonance parameter 
related to the link between the nodes i and j. To each node i attach a simply laced Lie 
algebra gj. Produce a reduced diagram Cji containing the node j by cutting the link adjacent 
to it in the direction i. Likewise produce a reduced diagram Cij containing the node i by 
cutting the link adjacent to it in the direction j. Then the Gc-theory decouples according 
to the rule 

= G (c-c tJ ) ® G^c-Cji)/ G(c-c tJ -c 3l ) ■ (3-3) 



(Tij — >00 



We depict this rule also graphically in terms of Dynkin diagrams: 

C o j-j — ► oo 



& Si 

C Cji C Cij / C Cij Cji 



gi gj 

According to the GKO-coset construction [Q] , this means that the Virasoro central charge 
flows as 

c gc ~^ c sc-c i:j c Sc-c ji ~ Cgc-c^-Cji ■ (3-4) 

The rule may be applied consecutively to each disconnected subgraph produced according 
to the decoupling rule ( |3.3| ). Note that this rule describes a decoupling and not a fusing, as it 
only predicts the flow in one direction and the limit is not reversible. From a physical point 
of view this is natural as also the RG flow is also irreversible. The rule (^^) generalizes 



a rule proposed in [26], which was based on the assumption that unstable particles are 



associated exclusively to positive roots of height two. 
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More familiar in the mathematical literature is a decoupling rule found by Dynkin [38] 
for the construction of semi-simple 2 subalgebras h from a given algebra g. For the more 
general diagrams which can be related to the g^-HSG models the generalized rule can be 
found in |3l|. These rules are all based on removing some of the nodes rather than links. 
For our physical situation at hand this corresponds to sending the masses of all stable 
particles which are associated to the algebra of a particular node to infinity. As in the 
decoupling rule ( |3.3| ) the number of stable particles remains preserved, it is evident that 
the two rules are inequivalent. Letting for instance the mass scale in gj go to infinity, the 
generalized (in the sense that gj can be different from At) rule of Kuniba is simply depicted 
as 

C nij — ► oo C — Cji C — Cjk 



Si Sj Sk Si Sk 

Clearly this can not be produced with (pT 



3.2 A simple example: The SU(4) 2 -HSG model 



(T23 the rule ( |3.3|) predicts the following flow 



We illustrate the working of the rule (3.3) with a simple example. We take g to be £J7(4), 
attach to each node simply an SU(2) algebra and to the links the resonance parameters 
o~i2, ci3i 0"23- This corresponds to the S'L r (4)2-HSG model. For the ordering 013 > o\2 > 



: = SU(4) 2 c = 2 

: = SU(3)f 2 /SU(2)2c = 1.9 
: = SU{3) 2 ® SU(2) V = 1.7 




Q'3 



• g = 5[/(2)f c=1.5 

«3 



The central charges are obtained from ( p. 28 ) using (3.4). Chosing instead the ordering 
o~23 > C13 > C12, we compute 

• • • g = SU(4) 2 c = 2 



oil a.2 03 

0-23 • • <g) • g = SU (3) 2 ®SU(2) X = 1.7 

Oil «2 «3 

<Ti3 is already decoupled 



^^12 • • (8i • g = 5(7(2)f 3 c=1.5 

ai «2 03 

It is important to note the non-commutative nature of the limiting procedures. For more 
complicated algebras it is essential to keep track of the labels on the nodes, since only in 
this way one can decide whether they cancel against the subgroup diagrams or not. 

2 The subalgebras constructed in this way are not necessarily maximal and regular. A guarantee for 
obtaining those, except in six special cases, is only given when one manipulates adequately the extended 
Dynkin diagram. 
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3.3 A non-trivial example: The (Eg)2-HSG model 

As by now we do not have a rigorous proof of the decoupling rule ([T^) , we take the support 
for its validity from the working of various examples. We will check below the analytic 
predictions of the rule against some alternative method. As the previous example was 
a simple pedagogical one, we will consider next a non-trivial one leading to an intricate 
prediction for the RG-flow. The confirmative double check below can hardy be accidental 
and we take that as very strong support for the validity of ( |3.3| ). 

We consider now the (Eg^-HSG model. In this case we have £ = 6, h = 12 such 
that we have 6 stable particles, 30 unstable particles and 15 resonance parameters. From 
the 5! possible orderings for the resonance parameters we present here only two concrete 
ones, which will predict different types of flows and mass degeneracies. Note that this 
degeneracy is not the unavoidable one resulting from the difference between the number of 
resonance parameters and non-simple positive roots that is 30 — 15, as discussed for ([ 



and (3.2). The degeneracies discussed here are a consequence of the particular choices of 
the resonance parameters. The conventions for the labeling of our particles are indicated 
in the following Dynkin diagram: 



a 2 



Ol\ 03 04 CK5 CVq 

We choose first the ordering and values for resonance parameters as 

0-13 = 100 > 0-34 = 80 > 045 = 60 > 0-56 = 40 > 0-24 = 20 . (3.5) 
According to the decoupling rule ( |3.3[) , we predict therefore the flow: 









Eq 


36 

T ~ 


5.14 


C16 


= 280 




SO(10)® 2 /SO(8) 


5 




Cl5 


= 240 




5O(10) ® SU(5)/SU(4) 


34 

T ~ 


4.86 


0-14 


C36 = 


180 


SO{8) <g> SU{5) <g> SU{3)/SU(4) ® SU(2) 


319 
70 


- 4.56 


0~12 


= 160 




is already decoupled 






0-35 


= 140 




SU{5) ® SU{4) ® SU(3)/SU(3) SU{2) 


61 

14 ~ 


4.36 


C26 


= 120 




SU{4)® 3 ® SU(3)/SU{3)® 2 ® SU{2) 


4.3 




0"13 


046 = 


100 


SU(4)® 2 ® SU(3) ® SU(2)/SU{3) ® SU(2) 


4 




0"25 


C34 = 


80 


SU{3)® 3 ® SU{2)® 2 /SU{2)® 2 


3.6 




0-32 


C45 = 


60 


SU{3)® 2 SU(2)® 2 


3.4 




<756 


= 40 




SU{3) ® SU{2) m 


3.2 




•724 


= 20 




SU{2)® 6 


3 





Note that eight particles are pairwise degenerate and we therefore expect to find 15 — 8/2 = 
11 plateaux in the flow. The first step which corresponds to one of these degeneracies occurs 
for instance at CT14 = 0-36 and we have to apply the decoupling rule twice at this point before 
we get a new fixed point theory. 
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Next we arrange the couplings as 

<7 45 = 100 > 0-34 = 80 > (7i3 = 60 > cr 56 = 40 > CT24 = 20 . 
and compute from fl3.3|) the flow 













5.14 


C16 


— ZoU 






5 




C15 


— Z4U 




C/Oflf^ QTTfK\ IQTT(A\ 


34 

- ~ 


4.86 


0"36 


= 220 




SU(5)® 2 ® SO(8)/SU(4)® 2 


33 

- ~ 


4.71 


0"35 


= 180 




SU{5)® 2 /SU(3) 


~ 4 51 

35 


0"26 


= 160 




SU{4)® 2 ® SU(5)/SU{3)® 2 


±££ ~ 4. 46 

35 


cr 14 


= <7 46 = 


140 


5C/(4)® 2 ® SU(3)® 2 /SU(2)® 2 ® S?7(3) 


4.2 




0"12 


= 0"25 = 


120 


5i7(4) 5^(3)® 3 /SC/(2)^ 3 


4.1 




0"45 


= 100 




S£/(4) <g> S£/(3)® 2 /SZ7(2) 


3.9 




0"34 


= 80 




S£/(3)® 3 


3.6 




C13 


= 0-32 = 


60 


5C/(3)® 2 (8) 5C/(2)^ 2 


3.4 




0"56 


= 40 




5i7(3) ® 5C/(2)® 4 


3.2 




0"24 


= 20 




SC/(2)® 6 


3 





(3.6) 



In this case we have only six particles pairwise degenerate and we expect to find 15 — 6/2 = 
12 plateaux. In the next section we find that the predictions made here are confirmed even 
for this involved case. 



4. How to detect unstable particles? 

In section 2 we described several arguments which predict the spectrum of unstable particles 
and now we will present some methods which allow to test these predictions. In particular 
with regard to the boostrap proposal this will be important, as it is not yet rigorously 
supported. Computing renormalization group (RG) flows will allow to detect the unstable 
particles. Roughly speaking, the central idea of an RG analysis is to probe different energy 
scales of a theory. We can flow from an energy regime so large that the unstable particle 
can energetically not exist to one in which it is formed. As a consequence, the particle 
content of the theory will be altered, which is visible in form of a typical staircase pattern 
of the RG scaling function. 

There are various ways to compute such scaling functions, such as the evaluation of 



the c-theorem [40] or an analysis by means of the thermodynamic Bethe ansatz (TBA) 
1 41]. In the first case we have to evaluate the expression 



c(r ) 



-1 I drr s (e( r )e(o)) 



(4.1) 



The main difficulty in this approach is to evaluate the two-point correlation function 
(O(r)O(O)) for the trace of the energy-momentum tensor depending on the radial dis- 
tance r. Most effectively, one can do this by expanding it in terms of form factors, for a 
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general recent introduction see e.g. [42] and references therein. It is well known that for 



many, even quite non-trivial, theories such form factor expansions converge extremely fast, 



see [27] for the compution of (fOj) for the SU(3)2-HSG model. 

Here we will concentrate more on the TBA, which is simpler to handle in most cases. 
As a prerequisite, one assumes to know all scattering matrices Sij(6) for the stable particles 
of the type i,j with masses m^rrij. Besides this dynamical interaction one also makes an 
assumption on the statistical interaction between the particles, which are choosen here 
to be of fermionic type. The TBA consists now of compactifying the space of this 1+1 
dimensional relativistic model into a circle of finite circumference R, such that all energies 
become discrete and functions of R. The function similar to (14-11), which scales now these 



energies takes on the form 

oo 

Ceff (r) = ^ m i J d0 cosh 9 m (! + e" £l(e ' r) ) . (4.2) 

« -oo 

One identifies the circumference R with the inverse temperature T and introduces the 
scaling parameter r = m/T, with m being an overall mass scale. The 6i(6,r) are the 
so-called the pseudo-energies which can be obtained as solutions of the thermodynamic 
Bethe ansatz equations 



cosh e = Ei{e, r) + Yypn * H 1 + e " £j )](#, r ) ■ (4-3) 



Here the * denotes the convolution of two functions (/ * g) (6) := l/(2n) J dO 1 f(9 — 8')g{6') 
and the S (for the stable particles only!) enter via their logarithmic derivatives ifij(0) = 
—idlnSij(6)/d9. The main difficulty in this approach is to solve ( |4.3j ), which are coupled 
non-linear integral equations and therefore not solvable in a systematic analytical way. 
Now it is clear, that the two functions (4.1) and (|4.2|) can not be the same, but 



nevertheless they contain the same qualitative information. The functions will flow through 
various fixed points, at which the theory become effectively conformal field theories and the 
normalizations are choosen in such a way that the values of both functions coincide with 
the corresponding Virasoro central charges. When the theory is not unitary, (|4.2| ) has to 
be corrected by an additive term to achieve this. Computing then a flow from the infrared 
to the ultraviolet, one passes now various CFT plateaux, where the changes are associated 
to the formation of unstable particles with mass ( |2.15 ). The challenge is of course to 



predict the positions, that is, the height and the on-set of the plateaux, as a function of the 
scaling parameter. The on-set is related to the energy scale of the unstable particles and 
thus simply determined by the formula ( |2.15 ). To predict the height is less trivial and the 



proposal made in [19] is that the decoupling rule ( |3.3| ) achieves this. It is important to note 
here that a — > oo in (|3.3| ), means in the RG context a 3> all other resonance parameters. 
In the follwing picture we present the numerical computation for the (Eg^-HSG model, 
which precisely confirms our analytical predictions made by the decoupling rule in section 
3.3 
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5,0 
4,5 

c(r) 

4,0 
3,5 
3,0 




a ,=\ OOjotf=80,n; =60,o =40,a„,=20 I CT i3 a 32 : 24 



"34 "«>"45 ""3 U 56 '"^24 



-160 



-140 



-120 



-100 -80 

log(r/2) 



-60 



-40 



-20 



Having confirmed the predictions of our decoupling rule with a TBA-analysis, let us 
now discuss how the results of this analysis are compatible with our bootstrap proposal with 
a simple example: We consider the processes ( 2.40|) , Q2.42 ). In order to be able to interpret 
the BW-formula for the production of the particle (123) from (12) + 3 or (23) + 1 one has to 
"make" (12) and (23) stable, which is achieved when o~\2 or 023 is zero. One has to do that as 
otherwise the BW can not be applied, it only makes sense for stable particles. The first not 
obvious result here is that the resulting mass for (123) turns out to be the same from both 
cases ( p. 42 ) (and in all other examples!!!). Looking at the outcome of the TBA calculation 
(see |19| ] for the numerics on this case) one finds precisely the value ( 2. 43] ) reproduced by 
the TBA at the onset ln(r/2) ~ — (T13/2 = —(a 12 + 023)/2. Now apparently in the TBA 
analysis a\2 or 023 are not zero, which seems to contradict the previous assumptions in 
the bootstrap analysis. To understand this, one should keep in mind the meaning of the 
steps in the TBA. The formation of the particle (123) takes place when its mass becomes 
greater than the energy scale of the RG-flow, i.e. when mexp(o"i3/2) > 2m/r. Let us chose 
for instance a 12 = 30, U23 = 60, then exp(o"i3/2) ~ exp(45) ~ 3.49 x 10 19 . To resolve the 
apparent contradiction, it is now important to note that the other unstable particles are 
formed several orders of mag nitude below at exp(30) ~ 1.06 xlO 13 andexp(15) ~ 3.72 xlO 6 . 
This means in comparison to the formation energy scale of particle (123) the parameters 
ci2> o"23 can be regarded as approximately zero, which is in agreement with the assumption 
in the bootstrap analysis! 

This is just the same picture as put forward in the decoupling rule: In the formulation 
we say 013 — > 00, but inside the TBA analysis this is a milder statement and just means 
(T13 S> 0"i2 ; 0"23- Further quite non-obvious confirmation comes from the results when 
choosing the parameters differently, i.e. in the example discussed here <Ti2 — > — <7i2- The 
two pictures completely coincide. With regard to previous studies, it is very important 
to note that the occurrence of the step at exp(o"i3/2) ~ 2/r had no explanation at all 
before. Only the onsets at exp(o"23/2) ~ 2/r and exp(o"i2/2) ~ 2/r could be explained as 
they correspond to the formation of unstable particles from two stable ones. The additional 



step (for other algebras there are far more) was a mystery pointed out first in [29]. In [19] we 
provided for the first time an explanation for this feature: We predict its height and on-set, 
thus explaining also why it is absent when the resonance parameters are chosen differently. 



-17- 



Integrable models with unstable particles 



For all other examples studied (not even all have been presented in this proceeding, see 



| IS] for more) this picture is completely consistent. 



5. Theories with an infinite amount of unstable particles 

We address now the question of how to enlarge a given finite particle spectrum of a theory to 
an infinite one. In general the bootstrap ( 2.1C| ), which is the central construction principle 



for the 5-matrix, is assumed to close after a finite number of steps, which means it involves 
a finite number of particles. However, from a physical as well as from a mathematical point 
of view, it appears to be natural to extend the construction in such a way that it would 
involve an infinite number of particles. The physical motivation for this are string theories, 
which admit an infinite particle spectrum. Mathematically the infinite bootstrap would 
be an analogy to infinite dimensional groups, in the sense that two entries of the S-matrix 
are combined into a third, which is again a member of the same infinite set. It appears 
to us that it is impossible to construct an infinite bootstrap system involving asymptotic 
states, although we do not know a rigorous proof of such a no-go theorem. Instead, we will 
demonstrate that it is possible to introduce an infinite number of unstable particles into 
the spectrum. 

5.1 q-deformed gamma functions and Jacobian elliptic functions 

In general, the .S-matrix amplitudes consist of (in)finite products of hyperbolic or /and 
gamma functions. Here we will argue, that to enlarge the spectrum to an infinite number 
one should replace these functions by q-deformed quantities or elliptic functions. Let us 
first recall some mathematical facts in this section. We start with some properties of q- 
deformed quantities, which have turned out to be very useful objects as they allow for 
instance to carry out elegantly (semi)-classical limits when the deformation parameter is 
associated to Planck's constant. Here we define a deformation parameter q and its Jacobian 
imaginary transformed version, i.e. r — * — 1/r, as 

q = exp(i7rr), q = exp(— iir/r), T = iK-\__i/Ki . (5-1) 

We introduced here the quarter periods K# of the Jacobian elliptic functions depending on 
the parameter t € [0,1], defined in the usual way through the complete elliptic integrals 
K t = jo W/2 (l - ^sin 2 e)- 1 ' 2 d6 . Then 

lim Ki = lim K\-e = tt/2, lim K\—> = lim Ke — ► oo . (5-2) 

It will turn out below that quantities in q will be most relevant for our purposes and 
therefore we state several identities directly in q, rather than q, even when they hold for 
generic values. The most basic q-deformed objects one defines are q-deformed integers 
(numbers), for which we take the convention 
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They have the obvious properties 



lim 



limfnl^ 

i-*0 H 

[n + mr]q 



n. 



1 

n/n' 



for m,m' ^ 
for m = in' = 



£->o [n' + rn,'r]g 

Next we define a q-deformed version of Euler's gamma function 

x ' [1 + "]|H4 



(5.4) 
(5.5) 



r^ + i) :=f] 



1 + n]q[n] 



(5.6) 



The crucial property of the function Tg, which coins also its name, is 

oo ^ / 1 -\- n\ x 
lim T s (x + 1) = lim IVx + 1) = TT — ^— ( ) = T(x + 1) . 

n=l ' 
We can relate deformations in q and q through 

g -(*+r/2-i/2)2 r q -(y)rg(l-j/) _ r g (-y/r)r g (l + y/r) 
g(»-H-/2-V2)a r 5 (a?)r 5 (l - x) Y q {-x/r)Y q {l + x/r) ' 

Frequently we have to shift the argument by integer values 

T q {x + l)=q x - l [x] q Tq{x) . 



(5.7) 



(5.8) 



(5.9) 



Relation ( |5.9| ) can be obtained directly from ( |5.6| ). As a consequence of this we also have 

m— 1 

Tg(x + m) = Tq(x) ^ x ' + ^ 1 [z + /] 4 mGZ (5.10) 

z=o 

m— 1 

r 5 (x) = r^(x-m) JJ ^~ 2 [x - Z - 1]^ meZ. (5.11) 



1=0 



Whereas ( |5.9| )-( |5TTo| ) hold for generic q, the following identities are only valid for q 



I>(l/2 - r/2)^(l/2 + r/2 
T 4 (x + 2r 



i=l 9 



r^(y + 2r 
rg(x i )r ( j(x i ± r/2 



(y<)r 4 ( W ± r/2 



hmfr r ^ ±r/2 



£ 1 /4 r „( 1 / 2 )2 

r«(y) 



(5.12) 
(5.13) 



ii m _L r *(— t -)r -(i - — ± - 

5-1 £1/4 « v 2^ 2 ; 9V 2iQ 2' 



i=i <?' 
1 

7T 



r^(xi) 



1=1 2=1 
P P 



for 



i=l 
x^O 



i=l 



- 19 - 



Integrable models with unstable particles 



Most of these properties can be checked directly by means of the defining relation (|5.6|). 
The singularity structure will be important for the physical applications. It follows from 
( |5.6| ) that the Tg-function has no zeros, but poles 



lim Tg(6 + 1) -> oo for m G Z, n G N . (5.16) 

0— >0p™=m-r-n 

Next we define 

r T ^ ,_ tanh(g - mx + g)/2 -j^- _ scfl, dnfl + 

-tanh(# + i7nr + ( 7)/2' ™ H W'-»*>M " sc + dn#_ (5 - 17) 

v " n=— oo ' 

with x G Q, cr G M and 6>± = (# ± Z7rx + a)iKi/ir. We employed here the Jacobian elliptic 
functions for which we use the common notation pq(z) with p,q G {s,c,d,n} (see e.g. [43] 



for standard properties). We derive important relations between the q-deformed gamma 
functions and the Jacobian elliptic sn-function 

sn(x) = i r«(^Ti)r«(i-ifc±a {5 18) 

i£* r,(i - 2^717^9(2^717) 

These relations can be used to obtain some of the above mentioned expressions. For 
instance, recalling that sn(i'Q) = 1, we obtain ( |5.12| ), With (5.6) we recover from this the 



well known identity sn(iKi_i/2) = i/l 1 ^. The trigonometric limits 

7T 

lim sn(x) = lim smx) = —7^—7 — = sin(x) (5.20) 

«-o v 1 <i^i r(f )r(i - f ) v ' y ' 

1 p(i 4. i^)r(i - fe) 

limsn(x) = lim sn(x) = - —. — —. — — = tanh(x). (5.21) 

i-+i { J g— 1 w i r(i-f)r(f) 

can be read off directly recalling ( |5.2[ ), ( |5.7| ) and presuming that ( |5.16| ) holds. We recall 
the zeros and poles of the Jacobian elliptic sn(#)-function, which in our conventions are 
located at 

zeros: = 2lK t + i2mK 1 _ e l,m€Z (5.22) 

poles: 9 l ™ p = 2lK t + i(2m + l)K\-i I, m 6 Z . (5.23) 

We have now assembled the main properties of the q-deformed functions which we shall 
use below. 

5.2 Generalizing diagonal S-matrices 



Here we follow [37] and propose a quite simple principle which introduces an infinite number 
of unstable particles into the spectrum. We note first, that in general many scattering 
matrices factorize in the following form 



S ab {9) = ST{0)S^(9). (5.24) 
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Here S^ n (9) denotes the so-called minimal S-matrix which satisfies the consistency rela- 



tions i)— vii) of section 2. The CDD-factor [18], only satisfies i)-vi) and has its poles in the 
sheet —it < Im.8 < 0, which is the "physical one" for resonance states. We note now that 
the minimal part is of the general form 



Sab(0) 



(5.25) 



with A being a finite set specific to each theory. Then we may define a new S"-matrix 

S ab (0)= J] Me,* ( 5 - 26 ) 



and note that the additional factor in ( 5.26| ) is just of CDD-type. Therefore ( 5.26| ) consti- 
tutes a solution to the consistency relations i)— vii) of section 2, and thus a strong candidate 
for a scattering matrix of a proper quantum field theory. Note that wheras ( 5.25| ) was a 
finite product of hyperbolic functions, the new proposal ( 5.26| ) contains, according to the 



identity ( 5.17 ) in addition elliptic functions, which lead to the desired spectrum of infinitely 
many unstable particles according to the principles outlined in section 2. 

5.3 Non-diagonal S-matrices 

We discuss now the elliptic sine-Gordon model, which may be related to the continuum 
limit of the eight-vertex model. The (anti)-soliton sector was studied many years ago in 



1 44]. In [45] we demonstrated that it is possible to associate a consistent breather sector 
to this model. Let us recall the argument by recalling the Zamolodchikov algebra for the 
soliton sector 

Z{0 1 )Z{0 2 ) = a(0 12 )Z(0 2 )Z(0 1 ) + d{0 l2 )Z{0 2 )Z{0 1 ) , (5.27) 
Z{0 l )Z(0 2 ) = b(0 l2 )Z{0 2 )Z{0 l ) + c(0i 2 )Z(0 2 )Z(0i) . (5.28) 

In comparison with the more extensively studied sine-Gordon model the difference is the 



occurrence of the amplitude d in (5.27), i.e. the possibility that two solitons change into 



two anti-solitons and vice versa. Invoking the consistency equations i)-v) one finds J44], 45 1 



a{d) = m n 



b{0) 

c{0) 
d{9) 



i>[ 



l+2fc 



A]I>[l-0 



l+2k 



A] 



l+2fc 



1±2* A] 



k=0 v r^-^A]i>[i + 6 -2 
rp[9-(k + i)\]r#>[i + e-k\] 



T e [-0 — (k + l)A]r« [1-9- k\] 



sn{i0/u) 



sn(i0/u + tt/p 
sn(7r/^) 



-a(0), 



a(9), 



sn{i9/v + tt/v) 
-Vlm(i6/v) m(v/u)a(6), 

r 4 [i + i]r { [-i]r 5 [i 



+ ± + $n[9 



re[i + e + z]r 4 [-9-$]T 9 [i + ± + r\r 4 [- 



(5.29) 



(5.30) 

(5.31) 
(5.32) 
(5.33) 
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Here we used A = —ti/K^v, 9 = i9/2K^u with v £ R being the coupling constant of the 
model. With regard to property vii), it is clear that it is important to analyse the singularity 
structure of the amplitudes ( 5,2S| )-( |5,32D to judge whether there exists a breather sector. 



For this we appeal to the relations ( 5.16j ), ( |5.22 ) and (5.23) and find the following pole 
structure inside the physical sheet 

9 n ™ p = 2mvK x _ l + i2nvK e , 9™ p = (2m + l)vK^ t + i(ir - 2nuK e ), 

9 l b ™ p = 2mvK 1 „ e + i(7T - 2lvK e ), 9 l b ™ p = (2m + l)uK X - t + i2lvK e , 

9 l ™ p = 2mvK 1 _ e + i2luK e , 9 l ™ p = 2mvK X -t + i(w - 2lvK e ), 

9 l £ p = (2m + \)vK X - t + i2l»K e , 9 l £ p = (2m + l)vK x ^ + i(ir - 2nuK e ). 

We took l,m £ Z,n £ N and associated always two sets of poles 9™™ p and 9™™ p to a(0), 
and 9fr™ p to 6(0) etc. One readily sees from this that if one restricts the parameter 
v > tt/2Ki all poles move out of the physical sheet into the non-physical one, where they 



can be interpreted in principle as unstable particles. This was already stated in [44|, where 
the choice v > tt/2K^ was made in order to avoid the occurrence of non-physical states. 
This is clear from our discussion of property vii) in section 2, as we would have poles in the 
physical sheet beyond the imaginary axis, which when interpreted with the Breit-Wigner 
formula leave the choice that either mg < or T 5 < 0, i.e. we either violate causality or we 
have Tachyons. The restriction on the parameters makes the model somewhat unattractive 
as this limitation eliminates the analogue of the entire breather sector which is present in 
the sine-Gordon model, such that also in the trigonometric limit one only obtains the 
soliton-antisoliton sector of that model, instead of a theory with a richer particle content. 
For this reason we relax here the restriction on u and note that the poles 

e bu P = 9 c2,P fOT 0<n<n max = [n/2uK e ],neN (5.34) 

are located on the imaginary axis inside the physical sheet and are therefore candidates 
for the analogue of the n tft -breather bound states in the sine-Gordon model. We indicate 
here the integer part of x by [x]. In other words, there are at most n max — 1 breathers for 
fixed v and I. The price one pays for the occurrence of these new particles in the elliptic 
sine-Gordon model is that one unavoidably also introduces additional Tachyons into the 
model as the poles always emerge in "strings". It remains to be established whether the 



poles ( 5.34 ) may really be associated to a breather type behaviour. 

Let us now see if the poles on the imaginary axis inside the physical sheet can be 
associated consistently with breathers. We proceed similarly as for the sine-Gordon model 
fiqj , even though in the latter approach the following ansatz is inspired by the classical 
theory and here we do not have a classical counterpart. We define the auxiliary state 

Z n (9 1 ,9 2 ) := ±= [Z(9 1 )Z(B 2 ) + (-l) n Z(9 l )Z(9 2 )} . (5.35) 

This state has properties of the classical sine-Gordon breather being chargeless and having 
parity (— l) n . Choosing thereafter the rapidities such that the state ( 5.35| ) is on-shell, we 
can speak of a breather bound state 

lim Z n (9 1 ,9 2 )= lim Z n (9 l ,9 2 ) = Z n (9) . (5.36) 



12 
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Here 9^ is the fusing angle related to the poles in the soliton-antisoliton scattering ampli- 
tudes. We compute now with the help of (5.27) and ( 5.2S{ ) the exchange relation 



z n (9 1 )z(9 2 ) = s bnS {e 12 )z(e 2 )z n (e 1 ) 



(5.37) 



where 



Mj7 -%+nK t ) 



S bn s(0) 



[£sn- 



, 7T 



-sn 



iO 7T 

— + — + nK t 
v Iv 



la 



(5.38) 



and 



_ = i> [i + 9 + \ - § ]i> \-e-\- f ]i> l-9 + ± + n jr r [f) 
r>[i-a ' A 



i + v + % - gji^-g - % - gjr^L-g + g + |JJ> 

, [1 _ g + A _ f]I> _ A _ | ]r> + A + n ]r> [_ 
-_l[0-H + A_^ A + Z] 2 2[ _g + n_A_ fcA _ Z] 2 
X $13^23 11 — ; ; " 

7 — 1 



4. A _ ni 
^4 2J 



•^+|-f] 



x 



X 



n I A _ 

2 ^ 



IT r / / A 

fc=0 i - * ~ 2 + 4 ~ 



oo 

1 I 1 z ; i 'H <- z 4 j_ 



■ 4 - jg- l - ■ z 4 - J<? 



Where $ 



'.J 



'.J 



with 0jj being the difference of the on-shell rapidities. What is remark- 
able here and can not be anticipated a priori, is that all off-diagonal terms vanish, thus 
as ( |5.37| ) expresses in the soliton breather scattering there is no backscattering. Similarly, 
but more lengthy, we compute the scattering amplitude between the n 4/i -breather and 
m^-breather 



Z n {9i)Z m {92) — S\ >n b m {9\2)Z m {92)Z n {9\ l 



(5.39) 



where 



s b n b m {0) 



1 - fen 2 -sn 2 ( — + (n + m)K e + - 



i0_ 
v 



(5.40) 



9 7T o / iO 



1 — £sn — sn ( h (n 4- ra)Kn 

v \ v 



si\{i9/u — tt/v + {n + m)K{) 
si\{i9/u + tt/v + [n + m)K() 
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and 



a = $13^14^23^24 



oo n—1 
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k=l 1=1 
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nn 
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oo m— 1 
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(5.41) 



The latter expression ( 5.41 ) is tailored to make contact to the expressions in the literature 
corresponding to the trigonometric limit. Also for this amplitude the backscattering is 
zero. 

The matrix S bn b m (9) also exhibits several types of poles, a) simple and double poles 
inside the physical sheet beyond the imaginary axis, b) double poles located on the imagi- 
nary axis, c) simple poles in the non-physical sheet and d one simple pole on the imaginary 
axis inside the physical sheet at 9 = 9 b = iv{n + m)Kg which is related to the fusing 
process of two breathers b n + b m — ► b n+m . To be really sure that this pole admits such 
an interpretation, we have to establish according to ( |2.11|) that the imaginary part of the 
residue is strictly positive, i.e. 



- i \im(9-9 b )S bnbm (9)>0. 
ti—>u b 



(5.42) 



The explicit computation shows that this is indeed the case, see [45]. 



Furthermore, It is very interesting to check if also (2.10) is satisfied for the fusing 
process b n + b m — ► b n+m . For consistency, all amplitudes have to satisfy the bootstrap 
equations 



s ib n + m (°) = s ib n {Q + ivmKi)S lbm (6 - ivnKi) , 



(5.43) 



for I G s, s} ;k,m + n < n max . Indeed, we verify with some algebra that ( 5.43 ) holds 
for the above amplitudes ( |5,38| ) and ( 5.40 ). 

Finally, we carry out various limits. Our formulation in terms of q-deformed quantities 
and elliptic functions is useful to make this task fairly easy. We state our results here only 
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schematically and refer the reader for details to [45|. We find 

1 / u^2nKi /ir+2imKi_£ / it 
/ 



elliptic sine-Gordon 



m ^ oj _> o 
I 



free theory 



T 

l/v —* ioo 



elliptic dJJL-ATFT 



m = 0,£ -» 
i 



sine-Gordon 



minimal -D^h^-ATFT 



\/v— *n 

Thus we can view the elliptic sine-Gordon model as a master theory for several other 
models. In the limit f->0we recover now all sectors, including the breathers, of the sine- 
Gordon model. The diagonal limit 1/u — ► 2nKi/ir + 2imK\_ijix is interesting as it yields 
a new type of theory, which we refer to as elliptic SO(2n + 2) = D^^-affine Toda field 
theory (ATFT). To coin this name for these theories seems natural as in the trigonometric 
limit we obtain from it the ordinary minimal -D^j^-ATFT. 



6. Conclusions 



We reviewed the general analytical scattering theory related to integrable quantum field 
theories in 1+1 space-time dimensions. We made a proposal for a construction principle of 
an S-matrix like object which describes the scattering between two unstable particles or an 
unstable particle and a stable one. We tested this proposal with various examples and found 
a remarkable agreement with the outcome of the thermodynamic Bethe ansatz in what 
concerns the particle content and the RG flow of the theories. We described the general 
Lie algebraic structure of theories with unstable particles and propose a decoupling rule 
which predicts the RG flow when some of the parameters in the theory become very large. 
Alternatively, we tested these analytical prediction with the TBA. Finally, we discussed 
how one can construct theories with and without backscattering which contain an infinite 
number of unstable particles. 
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